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Abstract 

In this paper, we investigate the existence and characterizations of the Frechet 
derivatives of the solution to time-harmonic elastic scattering problems with respect 
to the boundary of the obstacle. Our analysis is based on a technique - the factoriza- 
tion of the difference of the far-field pattern for two different scatterers - introduced 
by Kress and Pai'varinta to establish Frechet differentiability in acoustic scattering. 
For the Dirichlet boundary condition an alternative proof of a differentiability result 
due to Charalambopoulos is provided and new results are proven for the Neumann 
and impedance exterior boundary value problems. 

Keywords : Elastic scattering, Navier equation, Frechet derivative, far-field pattern, 
Dirichlet condition, Neumann condition, impedance condition, inverse scattering. 

1 Introduction 

The inverse obstacle scattering problem for time harmonic waves is to determine the shape 
of the boundary and the location of a scatterer from far field measurements of the total 
wave. This problem is of practical interest in some important fields of applied physics, 
as for example non destructive testing in linear elasticity. Although such an inverse 
problem is theoretically difficult to solve since it is ill-posed and nonlinear, one can apply 
numerical methods to recover an approximate solution. The use of regularized iterative 
methods via first order linearization requires the Frechet differentiability analysis of the 
far-field pattern of the solution to the forward problem with respect to the boundary of the 
scatterer. An explicit form of the first derivative is needed in view of its implementation 
in iterative algorithms. 

In acoustic scattering, Frechet differentiability with characterizations of the derivative 
as the far-field pattern of the radiating solution to a new exterior boundary value problem 
were investigated by Hettlich |10) and Kirsch |16| via variational methods, and by Hohage 
|12j via the implicit function theorem. These characterizations allow the numerical imple- 
mentation of the derivatives from the knowledge of the boundary values of the scattered 
wave only (see [51 [121 ES])- By the use of boundary integral equation methods, one can 
express the far field pattern of the solution to scattering problems in terms of products 
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of boundary integral operators with singular Schwarz kernels. The Frechet differentia- 
bility analysis of the far field was then developped by Potthast |25| [27] for the Dirichlet 
and Neumann acoustic problems via the Frechet differentiability analysis of the bound- 
ary integral operators involved, in the framework of Holder continuous and differentiable 
function spaces. A characterization of the derivative can be obtain by directly deriving 
the boundary values of the solution. This approach was extended to electromagnetism by 
Potthast [26J for the perfect conductor problem and to elasticity by Charalambopoulos 
[2] for the Dirichlet scattering problem only. More recently, the Frechet differentiability 
of the class of boundary integral operators with pseudohomogeneous hypersingular and 
weakly singular kernels - which includes the usual boundary integral operators occuring 
in time-harmonic potential theory - was analyzed by Costabel and Le Louer [ll[5l|22j, in 
the framework of Sobolev spaces. The analyticity of the integral operators with respect 
to the boundary is proven. As a consequence, it yields the possibility to establish the 
Frechet differentiability of the far-field pattern for any scattering problem by the use of 
boundary integral representations and the chain and product rules. In this way we obtain 
an additional implementable formula to compute the Frechet derivatives of the far field 
by deriving the boundary integral operators (see |14 | I15 | [19]). 

This paper is devoted to the Frechet differentiability analysis of the far-field pattern of 
the solution to elastic obstacle scattering problems in three-dimensional homogeneous and 
isotropic media via an alternative technique introduced by Kress and Pai'varinta in [18] 
to establish Frechet differentiability for sound-soft and sound-hard obstacles. It is based 
on repeated uses of Green's theorem and a factorization of the difference of the far-field 
pattern of the scattered wave for a fixed obstacle and a perturbed obstacle. An interesting 
feature of the method is that it only requires the continuous dependence of the boundary 
values of the solution on the boundary in order to prove the Frechet differentiabilty of 
the far-field pattern. This approach was extended to the perfect conductor problem by 
Kress in [17] and the impedance problem both in acoustic and electromagnetic scattering 
by Haddar and Kress in |7j. 

The paper is organized as follow : In section 2 we recall elementary results on time- 
harmonic Navier equations in Sobolev Spaces, following the notations of |T]. More details 
can be found in j6l [T^ . The far field identity for elastic waves in the case of a Dirichlet 
boundary condition was established by Alves and Kress in [1]. In section 3, we use this 
identity to give an alternative proof of the differentiability result due to Charalambopoulos 
and we improve the boundary condition satisfied by the Frechet derivative of the solution. 
In section 4 we apply this method - following ideas of Haddar and Kress - to establish the 
Frechet differentiability of the boundary to far field operator simultaneously for the cases 
of Neumann and impedance boundary conditions and again we provide a characterization 
of the derivative. 

2 The Navier equation 

The propagation of time-harmonic elastic waves in the three-dimensional isotropic and 
homogeneous elastic medium characterized by the positive Lame constants /i and A and 
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the density p is described by the Navier equation 

div a('u) + /jcj^ix = 0, (2.1) 
where cj > is the frequency. Here 

a{u) = A(divw)l3 + 2^z{u) and z{u) = ^([Vw] + [Vu]""") 

denote the stress tensor and the strain tensor respectively. Notice that I3 is the 3-by-3 
identity matrix and [Vu] is the matrix whose the j-th column is the gradient of the j-th 
component of u. We set A*u := div ct(u) = ^Aii + (A + /i) V div u. 

Let C be a bounded domain with a boundary T of class and outward unit 
normal vector n and let il"^ denote the exterior domain R^\0. We denote by H^{Q,), 
Hf^jyi'^) and H'^iT) the standard (local in the case of the exterior domain) complex 
valued, Hilbertian Sobolev space of order s G M defined on Q.'^ and T respectively (with 
the convention = L^.) Spaces of vector functions will be denoted by boldface letters, 
thus H' = {H''f. We set : 

H{n,A*) := {u £ H\n) : A*u £ L'^in)} , 

The space H{Q,A*) is an Hilbert space endowed with the natural graph norm. 
We use the following traces and tangential derivatives : 

d 

— — = n • V (normal derivative) , 
on 

d 

T(n,d) = 2fi— h Andiv+/in x curl (traction derivative), 

on 

d 

Ai(n,d) = — ndiv + n x curl (tangential Giinter's derivative). 

On 

The tangential gradient Vr and the surface divergence divp are defined for a scalar 
function u and a vector function v by the following equalities 



„ „ du ( dv 

Vn = Vrii + T— n, div v = divp v + n • -— 
On \ On 

and the tangential Giinter's derivative can be rewritten as follow: 

Mv = \^Yv]n - (divr v)n. (2.2) 

We note that, due to the trace lemma, U\y ^ H^{T) for it G H {0,, A*)LlH loci^"" , A*) . The 
normal derivative ^W|r the traction derivative Tu^j- are both defined as distributions 
in iJ~2(r) via the first Green formula (see [U [211 [23] lemma [271]). 

For two (3 X 3) matrices A and B whose columns are denoted by (01,02,03) and 
(61, b2,b3), respectively, we set ^ : i? = oi • 61 + 02 • 62 + " ^3- The following lemma is a 
consequence of the Gauss divergence theorem and the identity 

div (a(u) • v) = A*u ■ v + a{u) : e{v). (2.3) 
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Lemma 2.1 For vector functions u and v in H{Q, A*), it holds the first Green formula 

/ A*u ■vdx= I Tu-vds- I a(u) : z{v) dx. (2.4) 
Jo. Jr Jn 

The symmetry of the product o'(ii) : ^{v) = A(divti)-(divt)) + 2^£(ti) : e{v) = (y{v) : £(u) 

yields the second Green formula 

/ {u- A*v - A*u-v)dx = / {u-Tv -Tu-v)ds. (2.5) 
Jn Jr 

If u and V solve the Navier equation in then each term in (j2.5p vanish. 

Further we will use the following different representations of the traction operator T : 

Tu = o{u)-n = 2^M.u + {\ + 2^){(i\w u)n — X c\iy\u (2.6) 
f du \ 

= ( — + + (A + /Li)(div-u)n (2.7) 

du 

= {\ + 2fi)- AXw + (A + /i)n X curlit. (2.8) 

on 

Now we assume that the domain Vt has a connected boundary F. In the sequel we are 
concerned with the following exterior boundary value problems for elastic waves: Given 
vector densities / G H^{T) and g G H^^iV), find a solution u € Hioc{^'^,A*) to the 
Navier equation ()2.ip in il*^ which satisfies either a Dirichlet boundary condition 



u = f on F (2.9) 

or an impedance boundary condition 

Tu + iaujy/pu = g on F. (2-10) 

The impedance coefficient a is assumed to be a real non negative constant. The case 
a = yields the Neumann boundary condition. In addition the field u has to satisfy the 
Kupradze radiation condition 

lim r \ — iKpUp] = 0, lim r f — iKguA =0, r = Ixl, 

r^oo \ or I r— >oo \ or I 

uniformly in all directions. Here, the longitudinal wave is given by Up = — K~^Vdivu 
and the transversal wave is given hy Us = u — Up associated with the respective wave 
numbers Kp and Kg given by 

P~ X + 2p' '~ p ■ 

Solutions of the Navier equation satisfying the Kupradze radiation condition are called 
radiating solution. 
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The fundamental solution of the Navier equation is given by 

iKs\x-y\ / ^iKs\x-y\ ^iKp\x-y\ 

Irs + 





^ \ 4tt\x — y\ 4:1t\x - 

It is a 3 X 3 matrix and for j = 1, 2, 3 we denote by $j the j-th. column of ^. We have 
= ^{x,y)'^ = ^{y,x). From the second integral theorem ()2.5p . for a radiating 
solution u E Hiod^'^,^*) to the Navier equation (|2.ip . one can derive the Somigliana 
integral representation formula for x £ ^l'^: 

u{x) = {[TyHx, y)Vu{y) - $(x, y)Tyu{y)) ds{y), (2.11) 

where Ty = T{n{y),dy) and Ty^(x,y) is the tensor obtained by applying the traction 
operator Ty to each column of ^{x,y). For existence and uniqueness of a solution to the 
above boundary- value problems via boundary integral equation we refer to Kupradze [21J . 
The radiation condition implies that the solution has an asymptotic behavior of the form 

g2Kp I a:; I g^/^s |^| /" 1 \ 

U{x) = W^ix) H T-r-uf^x) + i-r-r) , |x| OO, 

X X \ / 

X 

uniformly in all directions x = -j — j-. The fields u'^ and it^ are defined on the unit sphere 

5^ in and known as the longitudinal and the transversal far-field pattern, respectively. 
We introduce the Z^-spaces 

Ll{S^) = {he L^{S^); h{x)-x = 0}, 
LliS^) = {heL\S% h{x)xx = 0}. 

We have uf G LI{S^) and tx^ e LjiS^). 



3 The exterior Dirichlet boundary value problem 

The scattering problem of time-harmonic waves by a bounded obstacle leads to special 
cases of the above boundary value problems. In this section we consider the rigid body 
problem. The total displacement field u + u^ is given by the superposition of the incident 
field w*, which we assume to be an entire solution of the Navier equation, and the scattered 
field u, which solves the Navier equation in Q^, the Dirichlet boundary condition 

u + It* = on r, 

and satisfies the Kupradze radiation condition. 

For x £ let W be the 3x3 matrix whose the j'-th column Wj is the radiating 
solution of (|2.ip and (|2.9p for the boundary value 

f = -'^j{x,-) onT, 
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and set V{x,-) = ^{x,-) + W{x,-), that is W and V are the scattered and the total 
field, respectively, for the scattering of a point source located at x G il.^. We note that V 
satisfies the reciprocity 

V{x, y) = \V{y, x)]^ , x, y € Q'', x 7^ y, 

which can be derived from the second Green formula ()2.5p . the Somigliana integral rep- 
resentation formula (|4.2p and the symmetry of the fundamental solution. 

Lemma 3.1 The unique radiating solution v £ Hiod^'^,^*) of the Navier equation (|2.ip 
satisfying the boundary condition (|2.9p for any f G i?2(r) admits the following integral 
representation 

v{x) = ^ [TyV{x, y)]^f{y)ds{y), x e Q''. (3.1) 

Proof. From the second Green formula (|2.5p on for the radiating solutions W and v 
we can write 

^ {[TyW{x,y)]'v{y) - [W{x,y)]''Tyv{y)) ds{y) = 0, 

for all X G Q^. Using the boundary condition for W and v, the symmetry of <I> and the 
Somigliana integral representation formula for v we obtain 

v{x) = ^ ([T,f (x, y)Vv{y) + [W{x, y)VTyv{y)) ds{y) 

= ^ ([r,ci>(x, y)]'v{y) + [TyWix, y)]'v{y)) ds{y) 

= ^ [TyV{x, y)]" v{y)ds{y) = ^ [TyV{x, y)V f{y)ds{y). 

■ 

For a fixed incident field w*, we consider the boundary to far field operator 

which maps the boundary of the rigid scatterer onto the far-field patterns u"^ and u"^ 
of the scattered field u. In order to describe the dependence of the operator F on the 
shape of the boundary F, we choose a fixed reference domain Q and we consider variations 
generated by transformations of the form 

X I— X + ^(x) 

of point X in the space M^, where ^ is a smooth vector function defined in the neighborhood 
of F. The functions 6 are assumed to be sufficiently small elements of the Banach space 
'^^(F,]R^) in order that (I + 0) is a diffeomorphism from F to 

Te = (l + e)T = {x + e{x); X G F}, 
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so that the surface Tg is a connected boundary of class of a domain Qg. The mapping 
F : 6 ^ F{Tq) is well defined in the neighborhood of the zero function in 't^'^{T,M.^). We 
then analyze the Frechet differentiability of at = 0. We want to prove the existence 
of a linear and continuous mapping J-"'(0) : 'ta'^iT) — )• Z^(5'^) x Lp{S'^) such that we have 
the following expansion in x Lp 

^(0-^(0) =^'(0)e + o(||e|h^2), when ^ 0. 

To this end, since we can interchange the differentiation with respect to the boundary 
and the passing to the limit \x\ — )• oo, we will establish the Frechet differentiability of the 
scattered field away from the boundary F. 

By ng we denote the exterior unit normal vector to Tg and, in what follows, we 
will distinguish the quantities related to the exterior Dirichlet scattering problem for the 
domain Qg through the subscript 9. We use the following identity established by Alves 
and Kress in [Ij pp. 13. 

Lemma 3.2 Assume that 0. C 0,g. Then 

ug{x) - u{x) = - [ [V{x,y)]'Ty{ue{y)+u\y))ds{y), (3.2) 

for all X S rig . 

Proof. From the second Green formula ()2.5p for $(x, •) and ti* we can write 

^ ([Ty<^{x,y)]'u\y) - ^{x,y)Tyu\y)) ds{y) = 0, 

for all X ^Vf^. Using the boundary condition of u we then obtain 

u{x) = - j^^{x,y)Ty{u{y) + u\y))ds{y), xGfl'. (3.3) 

From the second Green formula (j2.5p for the radiating solutions W and ug in we can 
write 

^ ([TyW{x, y)]'ug{y) - [W{x, y)]'TyUg{y)) ds{y) = 0, x G %. 

From the second Green formula (j2.5p on 0,g\0, for W and u*, the last equation and the 
boundary condition of ug we have 

[W{x,y)]''Ty{ug{y)+u'{y))dsiy) 

(^[W{x, y)]''Ty {ug{y) + u\y)) - [TyW{x, y)]'' {ug{y) + u\y))) ds{y) 

(3.4) 

[W{x,y)]''Tyu\y) - [TyW {x , y)]'' u\y)) ds{y) 



^ {[W{x,y)]''Tyu\y) - [r,M^(x, y)] V(y)) ds{y), 
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for X £ Using the boundary condition of W and u on F, the integral representation 
(j3.3p of u and the second Green formula for radiating solutions we obtain 



-tx(x) = - l_W{x,y)Ty{u{y) + u\y))ds{y) 

'[Wix,y)]''Tyu\y) + [TyW (x , y)]" uiy)) ds{y) 



[W{x,y)]'Tyu\y) - [TyW {x^y)]" u\y)) ds{y) 
[W{x,y)YTy{ue{y) + u\y))ds{y), 
for X G Jl^. From (|3.3p we can write for Uq 

ug{x) = - [ <l>{x,y)Ty{ue{y) + u\y))ds{y), x G (3.5) 



We obtain the identity p.2p by combining the last two equations. I 

Remark 3.3 As Kress and Pdivdrinta pointed out in 1 18] in the acoustic case, the lemma 
\3.S\ remains valid when the domain J7 is not strictly contained in Q.q and W can he extended 
as a solution to the Navier equation in the exterior ofVig. By theorem 5.7.1' in 123] pp. 
169, this can he assured ifV is analytic and does not differ too much from 0,. In this 
case the last equality in ()4.4p follows hy choosing an open domain D such that QCi^lg C D 
and then applying Green's integral theorem first in D\0, and then in D\0,g. 

Lemma 3.4 Assume that F is analytic. Then the following expansion holds 

ue-u = 1^ [TyV{;y)]''Bu{y)ds{y)+o{\\e\y2), (3.6) 

in H{G,A*) for all compact suhset G of and 6 sufficiently small, where 

Bu = -{e ■ n) \-{n X T{u + v})) x n + —^ — (n • T{u + u^))n \ . 

\fi A + 2/i J 

Proof. We use similar arguments as in the proof of theorem 3.1 in [18j for the analogous 
acoustic case. We denote by 5^ and Kg the integral operators on the boundary Tg with 
singular kernels 2^(x,y) and 2[Tri;^{x,y)] respectively. The fundamental solution ^ is 
pseudo- homogeneous of class —1. It can be shown that these operators are bounded from 
-ff~2 (Fg) to itself (see [6] and [24J, pp. 176). From (|3.5p and the jump relations, it can 
be deduced that the traction derivative of the total field (ug + it*) solves the boundary 
integral equation 

{l + K'g- ir^Sg) T{ug + u\^J = 2{Tu\^^ - irju\^J. 
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The operator Kg is not compact, therefore we use a regularization technique (see |21) ) to 
modify the integral equation as below 

{Id + Be)T{ue + wjrj = 'i^e{Tv}^Ve " ^^"frj 

where {Id + Bq) = He {1 + Kg — irjSg), Tig is a strongly integral operator and Bg is a 
weakly singular operator. The new integral operator I + Bg : H~2{Tg) — )• H~2{Tg) is 
a Fredholm operator of the second kind which is invertible with bounded inverse. We 
use the transformation Tg which maps a function Ug defined on Tg onto the function 
ug o {I -\- 9) defined on T. Reducing the analysis in [4J to the continuity and not the 
differentiability, we can prove that the boundary integral operators TgBgTg^ and rgTigT^^ 
depend continuously on the deformation 9 E "^^(F) and so does the inverse Tg{I+Bg)~^TQ^ 
from the Neumann series. Since the incident field is analytic on the boundary Tg, we then 
deduce that the total field satifies 

lire {T{ug + u%^) - T{u + ^ 0, \\9\\.^2 ^ (3.7) 

Since F is analytic, by theorem 5.7.1' in [23] pp. 169 the total fields {u + u^) and W can 
be extended as a solution to the Navier equation across F in the exterior of H ilg. We 
also have 

||t, {T{u + i.^)|rj - T{u + u')\\^-i^^^ ^ 0, \\9\y2 ^ 0. (3.8) 

By lemma [22] and using Taylor's formula together with (|3.7p - (|3.8p . the boundary condi- 
tion for V and the smoothness of ^(x, •) up to F for x away from the boundary, it follows 
that 

/ [V{;y)]''Ty{ug{y)+u\y))ds{y)= [ [V{;y)]''Ty{u{y) + u\y))ds{y) + o{\\9\y2), 

in H{G,A*) for all compact subset G of W^. From this, the first Green formula (|2.4p and 
the boundary condition for W we obtain 

ug-u = -j |([Vyj] : (y{u + u')')^_^^^- pu;'^[V]^{u + u')'^xdy + o{\\9\y2) 
where 

ni = {y e ng;y ^ n} U {y e n;y ^ Qg}, 

and x{y) = 1 if y G and y ^ 0, and x{y) = ~1 if y G i7 and y Qg. Any z £ 0,q can 
be represented of the form z = y + t9{y), with y G F and t > 0. We have : 

(a) x{z)dz = {9{y) ■ n{y))ds{y)dt + o(||e||^2), 

(b) {u + u'){y + t9{y)) = o{\\9\y2) and V{; y + t9{y)) = o{\\9\\^2), 

d 

(c) V(u + It*) = —{u + It*) • + o(||(9||<:^2) and. 
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(d) fori = 1,2, 3, 



d 



[VVj] : a{u + u') = a{Vj) : [V(it + u')] = TyVj ■ —{u + u') + oi\\9 



The expansion (a) for the volume form is a well-known result and we refer to |18| 13 |9] 
for a proof. Approximating the integral over by an integral over T we obtain that 

ue-u = -j^ [T,y(,y)]T^i!^tI^(y)(0(y) . n{y))ds{y) + omV^-)- 

To conclude, we express the normal derivative of {u + it*) in function of its traction 
derivative. Since u + li* = on F then the tangential Gunter's derivative M{u + u^)\r 
vanishes. Using ()2.7p . we obtain 



[n X T{u + u*)) X n = /I X '^{'^ + '"*)^ ''^j 
and using (j2.8p we obtain 

5(m + It* 



n-r(M + w*) = (A + 2^) n 



dn 



Theorem 3.5 Let T he analytic. Then the mapping T : <^2(r, _^ lI{S'^) x -Lp(S'^) is 
Frechet differentiable at 6 = with the Frechet derivative defined for ^ E "^^(TjIR^) by 

where v'^ is the far-field pattern of the solution to the Navier equation in Q'^ that 
satisfies the Kupradze radiation condition and the Dirichlet boundary condition 

= • n) ( -(n X T{u + n*)) x n + — (n • T{u + w*))n ) . 

Proof. This is a direct consequence of Lemmas 13.11 and 13.41 ■ 



Although theorem 13. 5l was proven under the assumption that T is analytic, we expect 
that, proceeding as in [18j, the result is also valid for boundaries. 



4 The exterior impedance boundary value problem 

Now we consider the scattering problem by a cavity are an absorbing obstacle Q. The 
total displacement field u + u^ is given by the superposition of the incident field u*, which 
we assume to be an entire solution of the Navier equation, and the scattered field u, 
which solves the Navier equation in Q'^, the impedance boundary condition 

T{u + u'') + iaoj^{u + u'') = on F, 
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and satisfies the Kupradze radiation condition. 

For X £ ^l^ let W be the 3x3 matrix whose the j'-th column Wj is the radiating 
solution of ()2.ip and ()2.9p for the boundary value 

g = -{T<i>j{x,-) + iauj^<^>jix,-)) on T, 
and set V{x, •) = $(x, •) + W{x, •). 

Lemma 4.1 The unique radiating solution v £ Hiod^^, ^*) of the Navier equation (|2.ip 
satisfying the boundary condition (j2.10p /or any g G i3"~2(r) admits the following integral 
representation 

v{x) = - j^[V{x,y)Yg{y)ds{y), x G 0^ (4.1) 

Proof. From the second Green formula (I2.5p on Q'^ for the radiating solutions W and t; 
we can write 

^ ([(Ty + iawVP)^(2;. y)Vv{y) - [W{x, y)]" {Ty + iau:^p)v{y)) ds{y) = 0, 

for all X £ VL'^. Using the boundary condition for W and v and the Somigliana integral 
representation formula for v we obtain 

v{x) = j^(^[{Ty + iauj^)^{x,y)Yv{y)-[^{x,y)]^{Ty + iauj^)v{y)^ds{y) 

= / (- [^^v + ia^VpWix, y)] ^v{y) - y)]^ {Ty + iaa;Vp)t^(y)) ds(y) 
= - ^ y)]'^(r, + iaa;V^)i;(y)ds(y) = - [V{x, y)]" g{y)ds{y). 

■ 

Here again we use a parametrization of the boundaries in order to investigate the Frechet 
differentiability of the boundary to far field operator 

F : r ^ tx- = {uf,u^) G Ll{S^) x L^S^) 

which maps the boundary of the Neumann or impedance obstacle $7 onto the far-field 
patterns and u"^ of the scattered field u. Thus we will consider instead the mapping 
T : 6 ^ F{rQ). The following lemma give a factorization of the difference of the scattered 
field for two neighboring impedance obstacles. 

Lemma 4.2 Assume that Q <Z Q.q. Then 

ue{x)-u{x)= f [TyV{x,y) + iau^V{x,y)]^{ue{y)+u\y))ds{y), (4.2) 
for all a; G Og. 
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Proof. From the second Green formula ()2.5p for ^{x, •) and we can write 

[{Ty + iau^)^{x, y)]^u'{y) - ^>(x, y){Ty + iauj^)u\y)^ ds{y) = 0, 
for all X G fi^. Using the boundary condition of u we then obtain 

u{x) = J^[(Ty + iaoj^Mx,y)]'{u{y) + u\y))ds{y), xGQ'. (4.3) 



From the second Green formula ()2.5p for the radiating solutions W and ug in i7g we can 
write for x £ ^Iq 

f {[(T^ + iaa;Vp)VF(x, - [W{x, y)]^ {Ty + iauj^)ueiy)] dsiy) = 0. 

From the second Green formula (|2.5p for W and it* in the last equation and the 

boundary condition of Ug we have 

[ [{Ty + iaoj^)Wix, y)V{ueiy) + u\y))ds{y) 
\{Ty + iaLo^p)W{x, y)]^{ue{y) + u\y)) 

-[W{x, y)Y{Ty + iau^){ueiy) + M^y))) dsiy) (4.4) 
[{Ty + iauj^)W{x, y)]'u\y) - [W{x, y)]" {Ty + iaoj^)u\y)^ ds{y) 

[{Ty + iau^)W{x, y)]^u\y) - [W{x, y)]^{Ty + iau^)u\y)) ds{y), 



re 



for X £ Using the boundary condition of W and li on F, the integral representation 
of u and the second Green formula for radiating solutions we obtain 



/r 



-w(x) = IJ{Ty + taoj^)W{x,y)]^{u{y) + u\y))ds{y) 

\{Ty + iauj^)W{x, y)Vu'{y) + [W{x, y)]'^ {Ty + iauj^)u{y)^ ds{y) 
[[{Ty + iau;^)W{x,y)]'u'{y) - [W {x , y)]" {Ty + iau ^)u' {y)^ ds{y) 
[{Ty + iau:^)W{x, y)\^ {ue{y) + u'{y))ds{y), 
for X G il^. From (|4.3p we can write for Uq 

ub{x) = - [{Ty + iau^)^{x,y)]^{uB{y) + u\y))ds{y), x G fi^. (4.5) 

JVf) 

We obtain the identity (j4.2p by combining the last two equations. ■ 
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The remark 13.31 is still availaible here, so that we expect that the identity (|4.2p remains 
valid when the domain Q is not strictly contained in Qq. 
We denote by Hr the mean curvature of T defined by 

= divr n. 

Lemma 4.3 Assume that F is analytic. Then the following expansion holds 

ug{x) - u{x) = - y)]Bu{y)ds{y) + o(| \e\ |<:^2(r,]R3)) , (4.6) 

in H{G,A*) for all compact subset G of Vt^ and 6 sufficiently small, where 
Bu = divr {{0 ■ n){aj(u + it*) — iaujy/pn ■ {ut + lij)"""}^ 

+ pu^e . n) ({I - -) {ut + u\) + (1 - -^-^) (it„ + <)nj (4.7) 
+ iau:y/p{e ■ n) ^[Vr(n + ii')]n + ^^^^ (divr(u + ii'))n -{u + u')H^ , 

with 

""'^""^ " ^ ( "^^^^ " " 'X?2^''") "3 + ^ ng ([Vrii] + [Vrn]^) Hs, (4.8) 

Ilg = I3 — n • n"*", Ut = {n X u) X n and u^ = u ■ n. 

Proof. We follow the same procedure as in the proof of lemma [3l4l We denote by 
5*51 and Kg the integral operators on the boundary Tq with singular kernels 2$(x,y) and 
2\Ty(^{x,y)]^ respectively. Assume that a > 0. From (I4.5p and the jump relations, it 
can be deduced that the restriction to T of the total field {ug + it*) solves the boundary 
integral equation 

{l-Ke- iaSe) {ue + u')\r, = 2u\^^. (4.9) 

The operator Kg is not compact. By regularization technique (see |21|). here again, we 
can modify the above equation in order to obtain an integral equation of the second kind 
which has to be solved for the unkown Uq + u^ in iJa (F). From this new equation and the 
convergence of integral operators as — )• 0, it can be deduce that the total field satifies 

||re(ne + it*)|r, -(u + it*)|r||^i^^^^O, \\e\\.g2 ^ Q (4.10) 
Then we can write that 

[{Ty + iauj^)V{-,y)\^{ue{y) + u\y))ds{y) 

[{Ty + iau^)V{;y)]'{u{y) + u\y))ds{y) + o(| |0| |c^2), 



13 



in H{G,A*) for all compact subset G of il.'^. Notice that the outer unit normal vector n 
to r can be extended in a continuously differentiable function, denoted again by n, on a 
tubular 

BT = {z = y + tn{y)- y G T, t G [-T; T]} 

for some sufficiently small T. By the first Green formula (j2.4p together with the following 
expansion for y G F (see [ZllH]) 

n{y) ■ {neiy) - n(y)) = - n(y) • Vr(0(y) • n(y)) + o(||/i||.^2) = oi\\h\\^^^) 

and the boundary condition for W we obtain 

U0-u = J ^i^div (^a{Vj){u + u'■) + ^auJ^/p{Vj ■ {u + Ui))nj ^ ^ xdy + o{\\6\y2) 
where 

and x(?/) = 1 if y G ile and y ^ Vt and x(2/) = ~1 if y G and y rj^i. Approximating 
the integral over VLq by an integral over T we obtain that 

Uq — U 

= div ((r{Vj){u + u'') + iauj^{Vj ■ {u + Ui))n) {9 ■ n)ds{y) + o{\\9\y2) 

For j = 1, 2, 3, we have 

div (o-(V,)(w + w^) = a(yj) : [V(m + it*)] - • (w + w*)) 

= [Vr(V,)] : (n3a(w + u')Us) - puj^{Vj ■ {u + u')) 

-iau^ (^-^V, ■ {u + u^) + [VrV,] : ((it* + wj) • n"^)^ 

and 

div ((F,- • {u + Wi))n) = + ^r) {Vj ■ {u + u,)) 
Collecting the two above equalities we obtain 

div (^a{Vj){u + It*) + iaujy/p [Vj ■ {u + tij))n^ 

= [Vr{V,)] : {U^aiu + u')U3 - iau^ {{ut + ul) ■ n^)) (4 

-puj^ {Vj ■ {u + ii')) + iauj^Vj ■ + Hr^ {u + it') 

To conclude we have to express n3CT(u+it*)n3 and ^(it+u*) in function of the tangential 
derivatives of u. First, we note that 

n3(x{u + u')U3 = Adiv(w + 'u*)n3 + /in3 ([Vrit] + [Vru]^) 
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Then we use the identity (|2.6p together with the boundary condition of u, which gives 

div(it + ■u*) = — divrftx + n*) — i-— ~ — ioJJ>n ■ (u + u^). 
A + 2/i A + 2/x 

The identity (j2.7p yields 

[ n X -^{u + It*) j X n = — [Vr(w + w*)]n — i—uj^ [n x {u + u^)) x n, 

end the identity (j2.8p yields 

^ ■ A Of 

Substituing all the above identity in ()4.1ip we obtain the characterization ()4.7p . ■ 

Theorem 4.4 Lei T be analytic. Then the mapping T : '^^{r,R^)^L'^^{S^)xLl{S'^) is 
Frechet diff'erentiable at = with the Frechet derivative defined for ^ G ^^(r,M^) by 



J-'(0)e = V. 



oo 



where is the far-field pattern of the solution to the Navier equation in that 
satisfies the Kupradze radiation condition and the impedance boundary condition 

(T + iaujy/p)v^ 

= divr (^{(, ■ n){CTj(it + u^) - iauj^n ■ {ut + u\)^]^ 

+ pu^ii . n) (^(1 - ^) {ut + u\) + (1 - {un + <)n) 

+ iauj^{i ■ n) (^Vt{u + u^)\n + y^^{^ivy{u + u'))n - {u + u')Hr ) , 



where the symmetric tensor o"j(m) is given by (j4 
Proof. This is a direct consequence of Lemmas 14.11 and 14.3 



When a = 0, the integral equation (|4.9p is not uniquely solvable but one can prove that 
the total field {ug + u^) is the unique solution of an hypersingular boundary integral 
equation. Via regularization method on can prove that the estimation (|4.10p is still valid 
in this case. For the Neumann problem, we then obtain the characterization: 



Theorem 4.5 Let T be analytic. Then the mapping T : <^^(r, R'^) ^ /.^(S"^) x ^^(S'^) is 
Frechet differentiate at 6 = with the Frechet derivative defined for ^ G ^^(r,M^) by 

where v"^ is the far-field pattern of the solution to the Navier equation in 17^ that 
satisfies the Kupradze radiation condition and the Neumann boundary condition 

Tv^ = divr ((^ • n)G^{u + it*)) + puj'^{i ■ n){u + u*), 

where 

^^u) = -^(divrit)n3 + A^n3 ([Vrtt] + [VvuV) Hg. (4.12) 
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